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The low-lying bound states of a microscopic quantum many-body system of n particles and the 
related physical observables can be worked out in a truncated n-particle Hilbert space. We present 
here a non-perturbative analysis of this problem which relies on a renormalisation concept and work 
out the link with perturbative approaches. 
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The ab initio construction of a rigorous quantum 
many-body theory able to describe bound particle sys- 
tems like molecules, atoms, aggregates and atomic nuclei 
has developed over a long period of time starting in the 
sixties [0. The quantities of interest are the spectro- 
scopic properties of low-lying many-body states of the 
systems, their energies and wavefunctions from which it 
is possible to obtain informations about relevant phys- 
ical quantities like multipole moments, electromagnetic 
transition rates, spectroscopic factors and other experi- 
mentally accessible observables. 

In practice the explicit resolution of the problem ne- 
cessitates the diagonalisation of the many-body Hamilto- 
nian in Hilbert space which is spanned by a complete set 
of basic states, in principle of infinite dimension, in any 
case generally very large. The information of interest is 
however restricted to the knowledge of a few energetically 
low-lying states which often possess collective properties. 
Hence it would be convenient to work in a finite truncated 
subspace of the original Hilbert space. Since the physical 
properties of the system should not depend on the dimen- 
sionality of the space of states the original Hamiltonian 
has to be redefined and gets an effective Hamiltonian 
which acts in the restricted space. The modification cor- 
responds to a renormalisation of the fc-body (generally 
k—1,2) interactions which govern the system. 

Many different renormalisation procedures have been 
proposed and successfully applied in the framework of the 
shell model and related microscopic descriptions. Rig- 
orous projection methods lead to effective Hamiltonians 
which can in principle be explicitly generated by means 
of perturbation techniques [Q, |[ ^, ^, ^ . Unfortunately 
there exists no straightforward control on the conver- 
gence properties of the perturbation expansions, espe- 
cially not when the interaction between the particles is 
strong as it is the case in atomic nuclei for instance 
Many attempts have been made in order to overcome this 
problem |0, ^, |^ . More pragmatic procedures have also 
been introduced such as two-body matrix elements cor- 
responding to truncated perturbation expansion contri- 
butions whose use may be justified by means of physical 
arguments llQ] . Even more phenomenologically the inter- 



action has been parametrised in terms of 2-body matrix 
elements which are fitted in such a way as to reproduce 
an arbitrarily selected number of experimental quanti- 
ties 10 . Very recently effective 2-body interactions have 
been constructed by means of a non-perturbative renor- 
malisation technique which cuts-off the large momentum 
components of the interaction Consequently this 

procedure eliminates the effects of the strong repulsion at 
short distances which affect essentially the high-excited 
states of the spectrum. 

All these methods are phenomenological and adapted 
to explicit calculations. However they do not tackle 
the problem related to the breakdown of perturbation 
methods in specific but frequently encountered situa- 
tions. This is the case in the framework of the projection 
method when two states belonging to different subspaces 
come arbitrarily close to each other [^ . 

The investigations which follow are developed in the 
spirit of former work based on renormalization concepts 
H, |l||6[ 0. We introduce a non-perturbative formu- 
lation of the bound state many-body problem. It grounds 
on the dimensional reduction of Hilbert space which gen- 
erates an effective Hamiltonian characterized by running 
coupling constants. The use of projection methods al- 
lows to link the present approach to perturbative meth- 
ods. We show explicitly the connection between the di- 
vergence of perturbation expansions and the concept of 
fixed points. 

Formal framework. Consider a system with a 
fixed but arbitrary number of bound particles in 
a Hilbert space H^^'^ of dimension N governed 
by a Hamiltonian H(^) ...gf where 

^gi^\ 92^\ ' ' ' J dp^^ 9^^^^ E^re a set of parameters 
(coupling constants) which enter the expression of 
the interaction operators in H^^\ The eigenvectors 
1*1^' (5'^-')) = 1, • • ■ , A^} span the Hilbert space and 
are the solutions of the Schroedinger equation 

i7(^)(g(^))|vl/f )(g(-'))) = A,(5W)|*f )(3(^))) . (1) 

The solutions are obtained by means of a diagonalisation 
which fixes both the eigenvalues \^Xi{g^^^), z = 1, • • • , A'} 
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and eigenvectors ||5'i^''(<7^^'')), i — 1, - ■■ , in terms 
of a linear combination of orthogonal basis states 
i = 1, . . . , N}. Since dim H^^^ = N is generally 
very large if not infinite and the information needed re- 
duces to a finite part of the spectrum it makes sense to 
try to restrict the space dimensions such that the phys- 
ical quantities related to the part of the system which 
is of interest remain the same as those obtained in the 
original space Ti^^K If the relevant quantities of interest 
are for instance M eigenvalues out of the set {K{g^^^)} 
then 

//(M)(,^(M))|^(Af)(^(Af))^ ^ \{g(M))\q,(^''\g(M))) (2) 



with the constraints 



(3) 



for i — 1,...,M. Eq. (3) implies relations between the 
sets of coupling constants g^*^-* and g^^"^ 
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(4) 



with k — l,...,p. The solution of this set of equations 
generates a new, effective Hamiltonian whose spectrum 
in the restricted space 7^*^*^^ is the same as the corre- 
sponding one in the original complete space. 

System at temperature T = 0. We now develop a gen- 
eral approach which allows to follow the evolution of the 
effective Hamiltonian of the system when the dimensions 
of the Hilbert space are systematically reduced. Using 
the Feshbach formalism |l^ we divide the Hilbert 
space n'-^^ into two subspaces, PH'-^^ and QH^^^ with 



dimPH*^) = - 1, dimg?^^"^ = 1 



(JV) _ 



(5) 



In the projected subspace PTt^^^ the system with energy 
E is described by the effective Hamiltonian 

H^ff{E) = PHP + PHQ{E-QHQy^QHP . (6) 



where we define 



H = Hn 



gHi 



(7) 



The set of basis states , i ^ 1, ■ ■ ■ , N} are the 

eigenvectors of Hq with the corresponding eigenvalues 
{ci, i = 1, • ■ • , N}. The expression Heff{E) is generally 
the starting point of theories which use perturbation ex- 
pansions [0. Here we proceed differently. We consider 



i=l 

which is the projection in PTi^^^ of an eigenvector 

N 

\^r) = j:-^\g^''')m 

i=l 



(8) 



(9) 



of Ti.^^\ If A^^^-* is the eigenvalue corresponding to l^'i^'') 
we look for the solution of 



(^)pl,T,W\ 



(10) 



We consider Pj^^^-*) to be the lowest eigenstate and 
QH^^^ to contain f. i. the highest one. We impose 
the lowest eigenvalue in the PH^^^ subspace to be the 
same as the one in the complete space 



(11) 



Projecting this expression on (^ij gives a relation which 
fixes g^''^~^\ the coupling constant which characterises 
the strength of the interaction in the PTi^^^ subspace in 
such a way that Eq. (0) is verified 



($i|i/.,,(AW)|P|*W) = Ar^(5W)a(^)(5(^)) . (12) 
Then 

($i|i/,/^(A(^))|P*W) = ^(}(^— )) . (13) 

The r.h.s. of Eq. (|l|) can be worked out explicitly 
by using given by Eq. (^) in the expression of 

H(.ff{E = aJ^''). The denominator in the second term 
of Heff{\\'^^) is a scalar quantity since dimQ'H^^) — 1 
where N designates the highest lying energy eigenstate 
of |$Ar). One gets 



^(}(A^— )) ^ Hf^^-)+HooA.(A(^)-7^^^) — ^(A|-oo) 

with 
and 

:^(^=($iii7(^-i)|p*(^') 



(14) 
(15) 



(16) 



and the transpose of the preceding matrix ele- 

ment. By construction the Hamiltonian depends 
on the new coupling constant g^^~^\ Imposing Eq. (|ll| ) 
leads to a second order equation for g^-^'^\ One gets 
explicitly 

where 

AN) 



N 



(18) 



(20) 
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and 



Fin — ^ a-i^ Hni 



(21) 



Since Eq. ( p7| ) is non-linear in g^^"^-* and has two so- 
lutions, g*-^^^^-* is chosen as the one closest to g*-^^ by 
continuity. The process can be iterated step by step by 
projection from the space of dimension TV — 1 to iV — 2 
and further, generating subsequently a succession of cou- 
pling constants g^*^^ at each iteration. At each step the 
projected wavefunction jP^^'^-') is obtained from l^^*'''). 
Going over to the continuum limit if N is large, one can 
derive a non-linear flow equation 

dg 1 ^dc ^ db ^ ^ d,a ^ ^^^^ 

dx 2a{x)g{x) + h{x) dx dx dx 

where a(a;), 6(x), c{x) and g{x) are the continuous exten- 
sions of the corresponding discrete quantities. Eq. ( p2[ ) is 
a non-linear differential equation which a priori can only 
be solved numerically. It has to be mentioned here that 
the resolution of the discrete or continuous flow equation 
requires the knowledge of the lowest eigenvalue and the 
corresponding exact eigenfunction in the k - dimensional 
spaces. The eigenvalue may be experimentally known. 
The eigenvector components are not directly accessible in 
practice. Information can be obtained from observables 
like the quadrupolc moment of the system in its ground 
state. They can be fixed rigorously , f.i. by means of a 
Lanczos diagonalisation procedure which determines the 
lowest eigenvalue and eigenvector. 

The present analysis shows how renormalisation con- 
cepts enter the problem of space reduction in the many- 
body problem. 

System at finite temperature T. We consider now the 
case of a quantum many-body system at temperature 
T = If the dimension of the Hilbert space H'--''^ is k 
its canonical partition function reads 



True 



(23) 



if the Hamiltonian contains p coupling constants. Going 
over to HC'^"^) 



one imposes 



(24) 



In the limit where [3 gets very large the trace reduces 
to the contribution of the lowest eigenvalue A^^*^"* |2^ 
constrained to stay constant through the dimensional re- 
duction process induced by Eq. (p4|). 
Using the Trotter formula 



(25) 



for n going to infinity and recalling Eq. (0) the trace 
defined in Eq. ( ^ ) gets a polynomial of order n in g^^'> 



(26) 



where i 



Eq. 



i ' 

I) leads to 



■ ,n| 



are fixed real coefficients. 



i=0 



i=0 



(k-l)gik- 



(27) 



Going to the continuum limit fc i— > x and taking the 
derivative of the partition function leads to the flow equa- 
tion 



dgjx) 
dx 



En 
1=0 



dhi [x) 



YTi=Qih,{x)g' ^{x) 



(28) 



Since n has to be very large in practice one may go over 
to an integral formulation of the discrete sums. 

Exceptional points and fixed points. As a last interest- 
ing point we finally address the problem concerning the 
divergence of perturbation expansions and the existence 
of fixed points of the running coupling constant g{x). 

It has been rigorously established that the eigenvalues 
\k{g) of H{g) — Hq + gHi are analytic functions of g 
with only algebraic singularities [20| . They get singular 
at so called exceptional points g ~ Qe which arc first order 
branch points in the complex g - plane. Branch points 
appear if two (or more) eigenvalues get degenerate. This 
can happen if g can take values such that Hkk = Hii 
which corresponds to a so-called level crossing. As a 
consequence, if a level belonging to the PTL subspace 
defined above crosses a level lying in the complementary 
QTi subspace the perturbation development constructed 
from Heff{E) diverges [|l^. Exceptional points are de- 
fined as the solutions of 



/(A(ge)) - det[H{g,) - X{g,)I] - 



and 



dX 



U = A(3e) 







(29) 



(30) 



where f{X{g)) is the secular determinant. It is now pos- 
sible to show that exceptional points are connected to 
fixed points corresponding to dg/dx = in specific cases. 
If {Xi{g)} are the set of eigenvalues the secular equation 
can be written as 



N 



n (A - A,) = 



(31) 



Consider A = Ap which satisfies Eq. (P9|). Then Eq. ( |3(i| ) 
can only be satisfied if there exists another eigenvalue 
Xq — Xp, hence if the spectrum is degenerate. This is the 
case at an exceptional point. Going back to the algorithm 
described above consider the case where the eigenvalue Ai 
gets degenerate with some other eigenvalue A^'^'' {g = ge) 
at some step k in the space reduction process. Since Ai 
is constrained to be constant, 



(32) 
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which is reahsed in any projected subspace of size k and 
I containing states |<i>i) and Going over to the con- 
tinuum Hmit and considering the subspaces of dimension 
X and X + dx one can write 



dx 







dXi{x) 
dx 



Consequently 



dXj dge 
dge dx 







(33) 



(34) 



Due to the Wigner - Neumann avoided crossing rule the 
degeneracy of eigenvalues is generally not fulfilled for real 
values of the coupling constant. There exist however spe- 
cific situations, like systems with special symmetry prop- 
erties |2^, ^ or infinite systems for which de- 
generacy can occur. In these cases Eq. (^) is realised 
if 



dge 
dx 



and 



dXj 

dge 



^0 



(35) 



The second relation works if crossing takes place and ge 
is a fixed point in the sense of renormalisation theory. 

Eq. ( |35| ) shows the connection between exceptional and 
fixed points in the framework of the present approach. 
Ground state degeneracy is indeed a signature for the ex- 
istence of phase transitions pT| , perturbation expansions 
break down at transition points. In fact, one identifies 
in the present result the properties of a quantum phase 
transition which is induced by level crossing ^3, p4| . 
The ground state wavefunction changes its properties 
when the (real) coupling constant g crosses the excep- 
tional point ge- There the cigenstates exchange the main 
components of their projection on the set of basis states 
\^,),i = l,...N. 

Co nclusions In summary, we developed a non- 



the physical constants which characterise them can be 
worked out. 
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